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ABSTRACT 

The kinematic evolution of axisymmetric magnetic fields in rotating 
magnetospheres of relativistic compact objects is analytically studied, based 
on relativistic Ohm's law in stationary axisymmetric geometry. By neglecting 
the poloidal flows of plasma in simplified magnetospheric models, we discuss 
self-excited dynamos due to the frame- dragging effect (originally pointed out 
by Khanna & Camenzind), and we propose alternative processes to generate 
axisymmetric magnetic fields against ohmic dissipation. The first process 
(which may be called induced excitation) is caused by the help of a background 
uniform magnetic field in addition to the dragging of inertial frames. It is 
shown that excited multipolar components of poloidal and azimuthal fields are 
sustained as stationary modes, and outgoing Poynting flux converges toward 
the rotation axis. The second one is self-excited dynamo through azimuthal 
convection current, which is found to be effective if plasma rotation becomes 
highly relativistic with a sharp gradient in the angular velocity. In this case no 
frame-dragging effect is needed, and the coupling between charge separation and 
plasma rotation becomes important. We discuss briefly the results in relation to 
active phenomena in the relativistic magnetospheres. 

Subject headings: black hole physics — magnetic fields — MHD — relativity 
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1. INTRODUCTION 

Interesting high-energetic phenomena have been observed in various compact 
astrophysical systems, such as pulsars, X-ray binaries and active galactic nuclei. Though 
the energy-release processes via radiation emission and jet formation have not been 
fully understood as yet, strong magnetic fields near the central objects can be a crucial 
component for explaining the observational features, and many works have been devoted 
to the developments of relativistic magnetohydrodynamical (MHD) models. In particular, 
Khanna & Camenzind (1994, 1996a) have recently proposed a self-excitation mechanism 
of axisymmetric magnetic fields, based on relativistic Ohm's law in Kerr geometry. This 
effect which is called self-excited gravitomagnetic dynamo is due to the coupling between 
relativistic frame-dragging of a rotating central object and rotational motion of surrounding 
plasma, and such dynamo action is expected to play an important role in the astrophysical 
phenomena as a trigger of relativistic plasma flows. Unfortunately, numerical calculations 
done by Brandenburg (1996) have shown that in a wide set of standard thin disk models 
around a rotating black hole no magnetic field can be maintained against ohmic dissipation. 
Cowling's antidynamo theorem for axisymmetric magnetic fields still holds even near a 
rotating black hole in the situations previously considered. Therefore, a new viewpoint will 
be necessary if the kinematic theory of resistive MHD presented by Khanna & Camenzind 
(1994, 1996a) is applied to the problem of generation of axisymmetric magnetic fields (see 
also Nunez 1997). 

In this paper we do not adhere to the accretion disk models, but we pursue analytically 
the processes permissible in rotating magnetospheres of compact objects. Though in the 
kinematic treatment of neglecting the feedback of magnetic fields due to the Lorentz force 
on plasma velocity the basic field equations become much simpler in comparison with the 
full MHD theory, some additional approximations are required to make any analytical 
approach possible. We would like to restrict attention to the simplified cases in which only 
the essential interactions between poloidal and azimuthal magnetic fields for allowing the 
existence of growing or stationary modes are preserved. 

The plasma injected into the rotating magnetosphere will partially accrete onto the 
central object and partially escape to infinity. The poloidal velocity of plasma flows, 
however, can remain sub-Alfvenic in some intermediate (quasi-equilibrium) region between 
the outer light cylinder and the surface of the central object, where the plasma angular 
velocity may be different from the Keplerian one, because gravitational forces do not 
dominate in comparison with any other interactions (see, e.g., Camenzind 1987; Takahashi 
et al. 1990). The stationary axisymmetric structure of the magnetosphere is mainly 
determined in the framework of ideal MHD theory under the frozen-in condition. We 
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consider the evolution of electromagnetic fields perturbed by the presence of small magnetic 
diffusivity. Though the motion of the plasma is assumed to be non-perturbed, a complicated 
dynamical evolution of electromagetic fields due to the poloidal flows will occur. Therefore, 
by fixing the poloidal velocity to be zero in the quasi-equilibrium region, we study the 
generation of magnetic fields which goes on slowly in a long diffusion timescale. Our 
purpose is to point out some basic aspects of the perturbed fields governed by the plasma 
rotation and the frame- dragging effect. We obtain the main results such that (i) if a 
background uniform magnetic field exists, it can sustain excited poloidal and azimuthal 
multipolar modes even in slow-rotation cases, and (ii) a sufficient charge separation 
generated by plasma rotation with a relativistic speed can cause self-excited dynamo 
without any frame-dragging effect. We discuss these processes of generation of magnetic 
fields (which have been missed in the above-mentioned numerical models) in relation to 
active phenomena observed in the relativistic magnetospheres. 

In the following we use units such that c = G = 1, and the axisymmetric stationary 
metric denoted by g a b has signs ( h + +). 



2. KINEMATIC EQUATIONS FOR AXISYMMETRIC DYNAMOS 

A contribution of accreting plasma and electromagnetic fields around a rotating 
compact object to the stationary axisymmetric gravitational field remains negligibly small. 
Therefore, we can always study the MHD interaction under a fixed gravitational field with 
the line element of the form 

ds 2 = g tt dt 2 + 2g t ^dtd<p + g^dcj) 2 + g„dr 2 + g ee de 2 , (1) 

where r, 9 and are the spherical coordinates, and the metric g a b is assumed to be 
independent of the time coordinate t and the azimuthal angle coordinate 0. The angular 
velocity of the dragging of inertial frames is denoted by uo = —gt^jg^- We define g to 
be the determinant of (g a b), and for the lapse function a = y—gtt + (g 2 ^/ fiV></>) we have 
\f— ~9 = ®-\/94>4>9rr9ee- We do not limit the metric to the Kerr form for later discussions. 
Further, we do not use explicitly the 3+1 formalism developed by Thorne, Price, & 
Macdonald (1986), and we treat relativistic Ohm's law in a covariant form written by 

F ab u b = 4nr](j a - Qu a ) , (2) 

where r\ is the magnetic diffusivity. According to the kinematic MHD theory the plasma 
4- velocity u a (r, 9) is also fixed, and Q = —j a u a is the electric charge density measured by an 
observer comoving with the plasma. The electric current density j a should be rewritten by 



-4- 



the electromagnetic field F ab via the Maxwell equations A-Kj a = V b F ab , where V& denotes 
the covariant derivative with respect to the metric g a b. 

Because we consider time-dependent fields under the assumption of axisymmetry, the 
poloidal magnetic components F r( p and Fq^ and the azimuthal electric component F t <f, are 
given by the single scalar potential r, 9) via the equation 

F a4> = d a ^ . (3) 

Then the four field variables F rt , F et and F r ' e remain to be solved, and the equation 
added to relativistic Ohm's law is the azimuthal part of the Faraday law 

d t F r0 + d r F et + d e F tr = . (4) 



These field equations for the kinematic evolution are still too complicated to discuss 
analytically the behaviors of solutions. Hence, the following investigation is limited to the 
models with no poloidal flow, i.e., 

u r = u 9 = , (5) 

which will be justified if plasma is located in a quasi-equilibrium region slightly distant from 
the surface of the central object, and a dynamical evolution of electromagnetic fields caused 
by the poloidal plasma flows with sub-Alfvenic velocities is not essential to the problem of 
dynamo action. Further, we assume f] to be a very small constant, and the time variation 
of fields is described by a long diffusion timescale such that t ~ r 2 /i]. Then it is convenient 
to introduce the variable T = r\t instead of t. As a result of this assumption concerning the 
order of rj the field equations can have consistent solutions, if the ratio of the amplitudes is 
understood to be 

rtFre/V = 0(1) . (6) 

Now let us give the field equations which are simplified according to the above- 
mentioned approximations. By virtue of equation (5) the poloidal part of equation (2) leads 
to 

F At u t + F A % = rjj A , (7) 
where the poloidal current density is approximately given by 

m A = -^e AB d B F , e re = -e dr = 1 , (8) 

in which the displacement current is neglected. (Hereafter the superscripts or subscripts 
written by A and B mean the coordinates r and 6.) Because we have Fa^ = Oa^, these 
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relations are used to express F At and F A ^ by ^ and F = r]y^gF r ' e , for example, in the 
approximated form of the proper charge density Q given by 



AttQ = --^=d b (^F ab ) ~ F^dAUt + F^dAU^, (9) 

which should be derived from equation (2) by using the current conservation V a j a = 
and the inequality \Q\ 3> \rjV a (Qu a )\. Note that for the stationary and axisymmetric 
velocity field u a we have d t u a = d^Ua = 0, while the rotational motion can generate the 
non-vanishing components &AUt and Oau^ for A = r,9. Then, except in the case that both 
u t and u<f, are constant, dbU a becomes non-symmetric under the permutation of a and b to 
assure the validity of equation (9) for the estimation of charge separation. 

From the Maxwell equations we also obtain approximately the azimuthal current 
density j^, in order to substitute it into the azimuthal part of equation (2) of the form 

u'dr^ = -MU - Qu4>) , (10) 

with equation (9) for the proper charge density Q. Then, we arrive at the final result of the 
evolution equation for \& 

d T m = S 1 + S 2 + S 3 , (11) 

where 



S 1 = _^d A (^— V*) , (12) 



— ^ d A ^(g 4>4> d A ^ + UtdAUcf, - u^AUt) , (13) 

e AB d A F{g H d B 0J - u^d B Ut + ud B u^)}. (14) 



(au^y^g 



The final term S3 can contribute to the excitation of ^ through the coupling to F. The 
ohmic diffusion is mainly due to the first term Si, and the role of S 2 (i.e., self-generation or 
self-destruction) will depend on the topology of ^. 

The Faraday law (4) is the evolution equation for F, which approximately reduces to 
the form 

d T F = ^{d A (^-d A F) - e AB d A *d B n} , (15) 



-g ' oP-vPg^ 

where Q = u^/u 1 is the plasma angular velocity. Note that the right-hand side of equation 
(15) is also decomposed into the terms for diffusion and amplification of F. In the following 
sections we will study the coupled equations (11) and (15) for \l/ and F to see the efficiency 
of excitation mechanisms. 



3. THE FRAME-DRAGGING EFFECT 



In the case of no poloidal velocity u A = we can give the specific energy and angular 
momentum of plasma denoted by —u t and as follows, 

-u t = 7{a 2 + g^uj(Q, - uj)} , = 73^(^-0;) (16) 

where 7 = — 1/ \Ja 2 — g^iVL — uj) 2 is the Lorentz factor of rotating plasma. If the plasma 
is co-rotating with the angular velocity of the background magnetosphere, the Lorentz 
factor 7 becomes very large in a region close to the light cylinder surface, and the term 
originated from the convection current density Qu^ will dominate in S3. In this section we 
would like to restrict attention to the frame-dragging effect, by neglecting any contribution 
of such charge separation under the condition Q = uj. 



3.1. The uj-Vt Dynamo 

Recall that in the numerical models of u$ = calculated by Brandenburg (1996) the 
dynamo action can work only if uj is taken to be artificially large. To see roughly this result 
from equations (11) and (15) with = 0, let us consider simplified forms of the metric 
components as functions of the coordinate 6 such that and — g are proportional to 
sin 2 ^, and uj, g rr and gee are independent of 6. (Now the Lorentz factor u f = 1/a depends 
only on r. This simplified metric may be regarded as the Kerr metric in slow-rotation 
approximation.) Then, by setting the time behavior of ^ and F to be exp(/iT), we obtain 

m {'- - LO* i)':A< -<),.F , (17) 



— x z a a 
and 

- L 2 )F = d 2 x F + aadj! , (18) 

1 — x z a 

where x = cos 9. The efficiency of dynamo action will be determined by the behavior of 
cr(r) dependent on the metric components as follows, 

sm9^/^g dr 



Further, L{§> and L 2 F in the left-hand sides of equations (17) and (18) represent the ohmic 
diffusion of ^ and F in radial direction, which are given by 



-7- 



and 



L 2 F = ^±d r ( * 9 8 r F) . (21) 

V^g ag<t>4>9rr 



We can decompose ^ and F into modes symmetric or antisymmetric with respect 
to the equatorial plane (Nunez 1996), and in this paper we would like to consider only 
the configuration of magnetic fields with the symmetry such that \&(— x) = ^f(x) and 
F(—x) = —F(x). (This corresponds to dipole-type topology of Quadrupole-type 
fields may be more easily excited near the equatorial plane. Here we do not pursue this 
possibility.) For physical modes the functions \& and F should satisfy the boundary conditon 
that both ^/(l — x 2 ) and F/(l — x 2 ) remain finite on the polar axis (x 2 = 1). Then, if the 
diffusion terms Li^ and L 2 F are neglected, it is easy to obtain 

q? = $ {l - (-l) n cos(ax)} , F = (-l) n a^f sm(ax) (22) 

as a stationary eigenmode with \i = 0. By virtue of the boundary condition for ^ and F 
the eigenvalue of a is given by a = mr, where n — 1, 2, • • •. Of course, if one takes account 
of a diffusion effect due to the terms and L 2 F, the minimum value of a should become 
larger than ir. However, for the Kerr metric on the equatorial plane, we can estimate the 
value of a to be 

2Ma(3r 2 + a 2 ) 
r 4 + a 2 (r 2 + 2Mr) " ' 1 } 

where M and a are the mass and rotation parameters, respectively. Therefore, one can 
expect no self-excitation of fields to occur near the Kerr black hole. In this sense the 
rotation of the black hole turns out to be too slow to excite the dynamo action. 



3.2. Induced Excitation 

Now let us propose an alternative process which can work even in the slow-rotation 
case and may be called induced excitation instead of self- excitation. The key assumption 
is the existence of a background poloidal field denoted by $b as a stationary solution 
of the vacuum Maxwell equations VbF ab = 0. (The typical example is given by Wald's 
(1974) solution for the Kerr hole immersed in a uniform magnetic field with the form 

— B {ag t( f ) + (#00/2) }.) If plasma is injected into the magnetosphere, the structure 
should be deformed by the motion of plasma. However, the original vacuum field in the 
background magnetosphere can remain dissipation-free and play a role of a stationary 
source field in equations (17) and (18). Then, we will be able to find a stationary solution 
written by ^ = ^> B + and F = F L . These perturbed parts * L and F L represent localized 
fields with amplitudes decreasing for large r, and the dynamical balance between the ohmic 
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dissipation and the excitation via the frame-dragging effect for ty L and F L is induced by 
the background field ^b- (This process has been also mentioned in Khanna & Camenzind 
(1996b) and Khanna (1997), and the numerical examples have been presented by Khanna 
(1998c).) 

To verify this induced excitation as a viable process, we write the stationary fields \I> 
and F satisfying equations (17) and (18) in the expansion forms 

* = (l-x 2 )tq2n(r) dP2n ; lix) , (24) 

dx 

and 

F = (1 -x ) 2^q2n+i(r) , (25) 

n=0 aX 

according to the boundary condition on the polar axis and the symmetry with respect to 
the equatorial plane. By the help of the recurrence relations for Legendre polynomials P n 
we have the equations for the coefficients q n as follows, 

9ee(Liq 2n ) ~ (2n + l)(2n + 2)q 2n = -(c 2 „+ig 2 n+i - Q2n-i52n-i) , (26) 

a 

and 

gee{L 2 q2n+i) - (2n + 2)(2n + 3)g 2 n+i = -acr(c 2n +2q2n+2 - c 2n q2n) , (27) 

where c n = (n + l)(n + 2)/(2n + 3). It is clear that the higher multipolar modes are 
generated from the lower ones through the action of the dragging of inertial frames denoted 
by o. 

Our main purpose is to point out a remarkable difference of the efficiency between 
the self-excitation and the induced one. Hence, for further analytical study, we consider a 
distant region where the metric components are approximately given by 

a = 1 , g H = r 2 sin 2 6>, g rr = g ee /r 2 = 1, (28) 

keeping the dragging of inertial frames written as 

uj = 2J/r 3 (29) 

for the angular momentum J of a central object. Equation (29) leads to 

a = 6J/r 2 < 1 . (30) 

In order to assure a = —r 2 du/dr to remain very small even for a small r in the following 
calculation, one may assume the behavior 

a = (r/r c ) 2 a c , a c = 6J/r 2 (31) 



- 9- 



in the inner region r < r c , where r c will be of the order of the radius of a central object. 

In the slow-rotation limit the recurrence relations (26) and (27) for n > 1 reduce to the 
approximated form 

An+l (n + 2)(n + 3) Un a °n n /ooA 

—^2 ~2 - y- 1 ) ' y 62 > 

because we can neglect q n+ 2 in comparison with q n in the right-hand sides. (Consistently 
with equation (32), the ratio q n +i/q n is assumed to be of the order of a, and the convergence 
of the expansion forms equations (24) and (25) is assured.) Now, if q n is known, it is easy 
to obtain the higher multipolar mode q n+ \. Note that in the flat spacetime the azimuthal 
component of magnetic field measured in orthonormal frame is given by B T = F/rjr sin 9. 
Then, q2n+i/r should vanish in the limit r — > oo and be regular in the limit r — > 0. The 
function g 2 n for any localized poloidal flux should satisfy the same boundary conditions. 
For a = we have the two independent solutions for each q n+1 as follows, 

q n+ i = r- ( " +2) , q n+1 = r n+3 , (33) 

violating either the outer boundary condition or the inner one. If there exists a localized 
field corresponding to a lower multipolar mode q n , however, the frame-dragging effect giving 
a 7^ can generate the higher one 

q n+ i = r-(" +2 ) [ r b n (p)p 2n+4 dp, (34) 

where 



o 



bn(r) = / {-l) n+l -^-Mp) d P ■ (35) 

Therefore, the key problem is the generation of the lowest mode qo, for which we obtain the 
equation 

d 2 q 2 6a 
ar r 5r z 

because c n = for n — — 1. Note that the remaining source field for q is only the lowest 
mode qi of the azimuthal magnetic field, satisfying the equation 

d 2 qi 6 2a 

d^T - ^ - ■ (37) 

Both modes should be self-consistently generated according to these coupled equations. 

As was previously mentioned, a localized solution as a result of self-excited dynamo 
will be prohibited for a small a. For example, we can check the suppression of the action 
even in the extreme case of assuming a sharp gradient in the angular velocity u> at r = r c : 
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The value of u is given by a nonzero constant Auo in the inner region r < r c , while it 
becomes zero in the outer region r > r c . In this case, though the localized modes q and 
qi can be continuous with values denoted by q c0 and q cl , the gradients dq /dr and dqi/dr 
must have discontinuous gaps estimated to be — 3g c o/r c and — 5g c i/r c , respectively. (We 
have q n ~ r n+2 for r < r c , while g n ~ r -0+i) f or r > rc _) Then, integrating equations (36) 
and (37) over the narrow region with the sharp gradients, we obtain the eigenvalue of the 
discontinuous gap given by r c Auj = 5\/5/2. 

Nunez (1997) has also treated a case with a sharp gradient in Q by using equation (29) 
for W, and he has claimed the existence of growing modes for a smaller discontinuous gap 
r c AQ < 1. This difference will be mainly due to the estimation of the diffusion term given 
by d 2l &/dr 2 , which was rewritten into the form e 2 d 2 ^/dx 2 under the scale transformation 
r = r c — (e/2) + ex. It seems to be unacceptable that the amplitude of the diffusion term 
is suppressed by the small factor e 2 . In fact, the straightforward application of the scale 
transformation will lead to d 2 "^ /dr 2 = e~ 2 d 2 ^ /dx 2 . In our calculation the diffusion terms 
d 2 qo/dr 2 and d 2 q\jdr 2 are estimated to be very large by virtue of the steep change of 
the gradients dqo/dr and dqi/dr. Then, the ohmic dissipation due to the diffusion terms 
dominate in the evolution, unless the discontinuous gap Aoj itself becomes unphysically 
large. 

Therefore, let us consider the induced excitation for the physically plausible behavior 
of a given by equations (30) and (31). In the case a = 0, the two independent solutions for 
go have the forms r 2 and r~ l corresponding to a uniform magnetic field and a dipole one, 
which will be relevant to the magnetospheres of black holes and neutron stars. Of course, 
these vacuum fields should be regarded as background fields in the magnetospheres, of 
which the origin cannot be attributed to the dynamo action discussed here. Our problem 
is rather to study the field generation against the ohmic dissipation in the background 
magnetospheres. We limit the analysis to the case of the uniform background field with the 
strength B . Then, for a ^ due to the frame-dragging effect the solution q is modified 
into the form 

% = ~y~ + v c Po ■ ( 38 ) 

To obtain q 1 in the first order of a c , it is sufficient to substitute the uniform background 
field into equation (37), and the lowest mode of localized azimuthal field is found to be 

* = * J #-i> (39) 

for r > r c , and 



c " c 



- 11 - 



for r < r c . This azimuthal part can work as a source of the localized poloidal field p in 
(36), and we obtain 

B rl^rl 45r c + 104 
675r r 3 4r 7 



ft = i^fG - + ^) (41) 



for r > r c , and 



£ or 2 4r 3 45r 4 



c w c 



for r < r c . The generated poloidal part po(r) has a maximum point in the outer region 
r > r c , and the poloidal field lines along which (1 — x 2 )p (r) is constant show a loop 
structure on the poloidal plane, which is maintained against the ohmic dissipation. In the 
slow-rotation limit such a modification of the background uniform field due to the generated 
poloidal field remains small. Nevertheless, we can expect that a remarkable structure of 
poloidal field lines as a result of the induced excitation appears in the magnetosphere, if the 
result presented here is extended to a fast-rotation case of the Kerr geometry. 

The role of the generated azimuthal field Bt = 3sin# cos 8qi(r) /rjr with the strength 
of the order of BoJ/r]r may be astrophysically more important even in the slow-rotation 
limit. In the presence of azimuthal magnetic fields, outflows of plasma from the central 
region may be efficiently drived by the Lorentz force. Further, a Poynting flux is excited 
if azimuthal magnetic fields exist in the magnatosphere, and in our approximation the 
Poynting flux vector P A (A = r, 9) is given by 

P A = ^F Bt F AB , (43) 

In the asymptotic region r ^> r c the components P r and P e are easily estimated to be 

pr = B^P 2 2 

27n?r 3 ' V ; 

and 

P e = --r^—r sin cos 0(1 + cos 2 #) . (45) 
47rr/r 4 

Interestingly, the Poynting flux propagates outward from the central region and converges 
toward the rotation axis along the curves given by (1 + cos 2 9)/r = const. Though the 
generation of azimuthal magnetic fields is also possible through ideal and non-relativistic 
MHD processes, the dynamical balance between the ohmic dissipation and the induced 
excitation due to the frame-dragging effect can be an origin of the magnetospheric structure 
responsible for producing high-energetic phenomena in the polar region. 
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4. 



THE EFFECT OF CHARGE SEPARATION 



In the case of vanishing no self-excitation of magnetic fields was found, even if an 
artificial sharp gradient of the frame- dragging angular velocity uo was assumed. This is 
mainly because a sufficiently large u is not allowed for any gravity around a central object 
rotating with a limited angular momentum. We also obtain the upper limit of the angular 
velocity f2 of plasma by the requirement (f2 — uo) 2 < a 2 / (see eq. [16]). Therefore, even 
in the cases Q ^ u, the u>-Q coupling will remain inefficient for self-excited dynamo, at 
least within the analytical framework developed here. (As was previously mentioned, the 
existence of growing modes claimed by Nunez (1997) may be a possible way to self-excited 
dynamo, if his estimation for the diffusion effect can be justified.) However, the situation 
may crucially change in the presence of charge separation given by equation (9). For the 
plasma co-rotating with the angular velocity of the background magnetosphere, the value 
of Uff, can become large without limit near the light cylinder surface, and a large gradient 
of u<j, (i. e., a large proper charge density Q) will be also allowed to occur there. Then, the 
azimuthal convection current Qu<f> which appears in equation (10) can play a key role for 
self-excitation of poloidal flux instead of the frame-dragging effect. 

To study the Q-Vt coupling as a mechanism of self-excited dynamo, let us restrict the 
following discussion to the case of no gravity and assume the plasma angular velocity to 
behave as Q = Q(R). (Hereafter, we use the cylindrical coordinates R, Z and 0.) Now the 
modes \l/ and F satisfying equations (11) and (15) can have the forms 



according to the assumed symmetry with respect to the equatorial plane. Using the Lorentz 



factor given by7 = l/yl — (RQ) 2 , the special relativistic versions of equations (11) and 
(15) reduce to 



= ip(R) cos(kZ)e' lT , F = f(R)sm(kZ)e flT , 



(46) 






(47) 



and 




(48) 



where L is the differential operator defined by 




(49) 



The plasma angular velocity fl may be nearly constant in an inner region (i.e., Q ~ fij, 
which is equal to the angular velocity of the background stationary magnetosphere with 
a rigid rotation), while it should decrease as R — > 1/f^. For mathematical simplicity, we 
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represent such a behavior as a sharp gradient in Q at R = R c < (Nunez (1997) has 

also studied this case in terms of the u-Q coupling without considering charge separation.) 
By virtue of the discontinuous gap AQ the gradients dip/dR and df/dR should also have 
the discontinuous gaps A(d^/dR) and A(df/dR) at R = R c . Then, using the equality 

HQ t 7/ j— \ i — \ 1 , 1 _ R£l s 

d 7 = -d(RQ + -In -) , (50) 

the integrations of equations (47) and (48) over the narrow region with the sharp gradient 
in f2 lead to 

^ , kmm + y_^ h (51) 

and 

A{-§-) = -tyAtm , (52) 
where ip = ip c and / = f c at R = R c . 

For further analysis we consider the discontinuous change from Q = Qi (i.e., 
7 = 1/^1 - RlVtj = 7i) to n = (i.e., 7 = 1). Then, for the modes with k > 1/R C , the 
amplitudes of ^ and / should decrease with the forms exp(—k (R — Rc)) at R > R c and 
exp(ki(R — R c )) at R < R c , where 



k Q = VA; 2 + /i, (53) 

and 

k t = y/P+^/i. (54) 

From these behaviors of ip and / we can easily estimate the discontinuous gaps A(dip/dR) 
and A(df/dR) at R = R C) and the growing rate /i is found to be 



H a(7i 



A; 2 6(7.) ' 

if the eigenvalue problem of equations (51) and (52) is solved, where 



(55) 



a = ^^H ll + f^i)-2, 



(56) 



and 

6 = 1 + 77TT (7i + v ^- + Tln(7 l + i v / 7 2 -l)-2 v /7,-l )1/2 ' (5?) 
The stationary mode with u = corresponds to the case 7$ = 7 m satisfying a(7 m ) = 0, and 
the numerical value is estimated to be 7 m ~ 5.55. 
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If plasma rotates with a relativistic velocity giving 7 > 7 m in a region around a central 
object, a sufficient decrease of Q within a range R c (l — 5) < R < i? c (l + 5) (5 <C 1) can 
excite growing magnetic fields with the scale of 1/k much smaller than the distance R c . 
An interesting possibility is that this self-excitation of magnetic fields produces relativistic 
outflows of plasma across the light cylinder as a result of the back-reaction, though it is 
a process beyond the scope of the kinematic theory considered in this paper. The spatial 
variation Ail of the angular velocity also may become smaller as the back-reaction works. 
Then, the self-excited dynamo will stop, and the acceleration of outflows will occur only 
through a stationary MHD process. The very active phase of self-excitation of magnetic 
fields and violent plasma acceleration (which will be responsible for flare-like events of 
radiation flux) can reopen only when a region with highly relativistic angular velocity 
appears again. 



5. CONCLUSIONS 

In a very simplified model we have succeeded in showing the self-excited dynamo 
through the coupling between fl and Q without any frame-dragging effect: Poloidal 
magnetic field generates azimuthal magnetic field by the help of differential rotation of 
plasma, and the azimuthal field induces charge separation if the Lorentz factor of the 
plasma rotation is not constant. The azimuthal convection current carried with the rotating 
charged plasma can amplify the original poloidal field. Though in our calculation a sharp 
gradient of Q has been assumed, the condition essential to the self-excited dynamo would 
be the existence of large spatial variations A(Ril) and A7 such that A(Ril)A'y > 4.5 
within a radial scale AR smaller than R. Different from the angular velocity uo of the 
dragging of inertial frames, we can expect A7 > 1 as an astrophysically permissible case, 
e.g., if we consider a rotational motion of plasma near the light cylinder in magnetospheres 
of relativistic compact objects. However, a siginificant poloidal motion may occur near 
the light cylinder, before 7 of plasma rotation becomes much larger than unity. Then, 
the charge separation may be suppressed in the magnetosphere. The self-excited dynamo 
due to the fl-Q coupling will be able to work only if strong background fields balance the 
centrifugal force in a plasma rotating with highly relativistic speed. Hence, this dynamo 
action should be understood to be an origin of flare-like instability permissible near the 
light cylinder rather than a mechanism of generation of strong background fields. 

We have also discussed the induced excitation of magnetic fields which occurs through 
the frame-dragging effect acting on uniform background magnetic field. This is a process 
leading to a new equilibrium configuration of magnetic field lines against the ohmic 
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dissipation. Our analysis has been limited to the case of slow rotation, in which the 
higher poloidal multipoles remain very small, and the extension to rapid rotation is an 
important problem to be solved. Though no flare-like activity is expected for this process, 
the magnetospheric structure with outgoing Pynting flux convergent toward the rotation 
axis is an interesting result. 

In this paper we have clarified only the fundamental aspects of kinematic generation 
of magnetic fields based on relativistic Ohm's law, and the astrophysical impications for 
black-hole or neutron-star magnetospheres should be confirmed in more realistic models. 
In particular, taking account of the presence of poloidal plasma flows will be an important 
task, which may work for antidynamo. Further, in the cases Q ^ u, the Q-Q coupling 
also can be an important origin of the induced excitation even if the self-excited dynamo 
does not occur. The combined effects of the u-Q and Q-Q couplings remain unclear in 
this paper. To treat more appropriately the problem of charge separation, one may need 
two-component plasma theories (see Khanna 1998a, 1998b). Axisymmetric dynamo action 
in charged plasma is an interesting problem in future investigations. 
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